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Abstract 

An analytic proof is given of the necessity of the Borland-Dennis conditions 
for 3-representability of a one particle density matrix with rank 6. This may 
shed some light on Klyachko's recent use of Schubert calculus to find general 
conditions for A^-representability. 



1 Introduction 

The recent announcement by A. Klyachko [8j of the solution of the pure state A^- 
representabihty problem for fermionic one-particle density matrix observes that this is 
the first new result since the work of Borland and Dennis [2\ in the early 1970's. There 
may therefore be some historical value in unpublished work of the author from that 
time, which makes a connection between the Borland-Dennis conditions and Weyl's 
problem. The latter asks for conditions on sequences {ak} , {bk} , {ck} which ensure 
that there exist self-adjoint matrices A,B,C with eigenvalues ak,bk,Ck respectively 
such that A + B = C. The first complete solution to Weyl's problem was given by 
Klyachko [7] in 1998. 

Let 7 be a density matrix normalized so that Tr7 = A^. The pure state 
A^-representability problem for fermions asks for necessary and sufficient conditions on 
7 for the existence of an antisymmetric A^-particle state whose one-particle reduced 
density matrix is 7. Let R denote the rank of 7. For the case A^ = 3 and R = 6, 
Borland and Dennis gave a pair of conditions on the eigenvalues of 7 which can 
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be written as follows under the assumption that they are arranged in non-increasing 
order. 

Ai + Ae = 1, A2 + As = 1, A3 + A4 = 1 (1) 

Ai + A2 < A3 + 1 (2) 

Note that ([T]) can be written compactly as A^ + Ay.^ = 1 for k = 1, 2, 3 

Borland and Dennis [2J proposed their conditions on the basis of numerical stud- 
ies and gave a proof of ([2]) under an assumption, which is equivalent to ([1]), about 
the pre-image of 7. In this note, we show that ([T]) is a necessary condition for 
iV-representability when = 3 and R = 6, completing the analytic proof of Bor- 
land and Dennis. We begin with some background material in Section [21 In Section [3] 
we present a proof of the necessity of ([1]). In Section H] we give a different, independent 
proof of the necessity of the inequality ([2]) from Weyl's inequalities. For completeness, 
we include a proof of sufficientcy of ([1]) and ([2]) in Section [51 In Sections [6] and [7] 
we present some partial results for the cases = 3 and i? = + 3 in the hope of 
providing some intuition behind the success of Klyachko's approach to a full solution. 



2 Notation and background 

In this note, we write the eigenvectors of 7 as \(/)k) so that 

1 = ^^k\<t>k){(f>k\- (3) 
k 

We will let A denote the anti-symmetrization operator and use the notation [fj, fk, fe] = 
Afj{xi)fk{x2)fe{x3) to denote a Slater determinant. The notation ( , )m indicates a 
partial inner product on a tensor product of Hilbert spaces. 

We need some results from Section 10 of Coleman's fundamental paper [3j. The 
first is Theorem 10.6 in [3] 

Lemma 1. (Coleman) The one-particle density matrix^ is N-representable with pre- 
image |\E') = V^i A |0i) ® |$i) + Vl ~ '^i |'^'2) if ctnd only if it can be written in the 
form 

7 = Ai|0i)(0i| + Ai7i + (1-Ai)72 (4) 

where 71 is the {N — l)-representable reduced density matrix of and 72 is N- 

representable with pre-image $2 satisfying 

(01, $2)1 = ($1, $2)2,3...iV = 0. (5) 
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The next two results are Theorems 10.2 and 10.4 respectively in [3j. (See also |10j.) 

Theorem 2. A one-particle density matrix 7 is 2-representable if and only if all non- 
zero eigenvalues are doubly degenerate. If there are no other degeneracies and the 
eigenvalues are written in non-increasing order so that X2k-i = > X2k+i, then the 
pre-image of 7 must have the form 

k 

Theorem 3. When N = 2n-\-l is odd and the one-particle density matrix j has rank 
R = N + 2, it is N-representable if and only if Xi = 1 and the remaining eigenvalues 
are doubly degenerate. 

3 Necessity of the condition + Xj-k = 1- 

To show that ([1]) is a necessary condition for 3-representability when R = 6, observe 
that since Ai = (0i,70i) it follows from (jlj) that 

(01,71 0l) = (01,72 0l) = 0. 

Therefore, 71 and 72 have rank < 5. It then follows from Theorem [3] that one can 
write 

72 = \9i){9i\ + \a\'^\92){92\ + |an^3)(^3| + \b\'^\9A){9A\ + \b\'^\95){95\ 

with |ap + |6p = 1 and [$2) = ci[9i, 92, 93] + b[9i, 94, 95]. There is no loss of generality 
in writing $1 = Ej<fc fi'fc]- 

We first consider the case in which both a,b ^ 0. Then a simple computation 
shows that ([5]) implies 

I'^'l) = C24[fi'2, 94] + C25[fl'2, 9^] + C34[fi'3, fi'J + CasfS'S, 9^] 

SO that (5^1, $1)1 = 0. Defining l^g) = \9i), gives Ae = 1 — Ai and one can rewrite (jll) 

as 

7 = Ai|(/)i)(0i| + (1 - Ai)|06)(06| + Ai7i + (1 - Ai)72 (7) 

where 72 = 72 — is the reduced density matrix of \Gi) = {gi, $2)3 = (i[92,93] + 

b[g4,g5]. Thus, in the orthonormal basis {92, 93, 94, 95} we find 

|C24P + |C25P C24C34 + C25C35 \ 

C24C34 + C25C35 |c34p+|c35p 

|C24P+|C34P C24C25 + C34C35 

\ C24C25 + C34C35 I C25 p + 1 C35 p y 
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The key point is that 71 is block diagonal and can be diagonalized by a block diagonal 
unitary transformation which mixes only within pairs §2,93 and 5^4, 5^5 leaving the 
Slater determinants in Gi unaffected (except possibly for a phase factor which can 
be absorbed into the new basis). Denoting the new basis as (pk, we now have IGi) = 
(^[02) (ps] + b[(p4, (^5]. Then either by explicit computation or from Coleman's proof [1] 
of Theorem [2], one can write = s[02, ^4] + ^[03, ^5] with |sp + |tp = 1. Thus, the 
eigenvalues of 7 satisfy 

A2 = Ai|a|2 + (1 - Ai)|s|2 (8a) 
A3 = Ai|6|2 + (1-Ai)|s|2 (8b) 
A4 = Ai|ap + (1 - Ai)|tp (8c) 
As = Ai|6|2 + (1-Ai)|t|2 (8d) 



which implies 



A2 + As = A3 + A4 = Ai + (1 - Ai) = 1. (9) 



We now consider the possibility that one of a,b is zero, in which case, [$2) is a 
single Slater determinant and there is no loss of generality in writing as $2 = [91,92, 93]- 
Then (|5]) implies that one can write 

I'^'i) = Yl Y ^jk[9j, 9k] +494,95] (10) 

j=l,2,3A:=4,5 

Now regard Xjk as a 3 x 2 matrix and observe when U, V are 3x3 and 2x2 unitary 
matrices, Y = UXV'^ corresponds to a basis change which mixes 91,92,93 and 9^,95 
among themselves. By the singular value decomposition we can find U, V such that 
only 1/24 and 2/35 are non-zero. Thus, in the new basis which we call (pk 

1*1) = y24[02, (pi] + 2/35[03, 05] + c[04, ^s] • (H) 

Again writing (p^ = 91, we find that the pre- image of 7 has the form 

1^) = ai23[(pl, (p2, 03] + a246[02, 04, M + a356[03, 05, M + a456[04, 05, M (12) 

which implies ([T]). 



4 Necessity of the inequality (0) 

We now prove that the inequality ([2]) is necessary for A^-representability. When 7 has 
the form ([3]) and ([1]) holds, its pre-image can be written in the form 

I ^) = a;ooo [01 , 02 , 03] + a;ooi [0i , 02 , 04] + a^oio [0i , 05 , 03] + a^oii [01 , 05 , 04] 

+ a;ioo[06,02,03] +a:ioi[06,02,04] + a;iio[06, 05, 03] + a:ill[06,0S,04]- (13) 
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In this form, there is no loss of generahty in assuming that the are arranged in 
non-increasing order. If we now define 



g I -^UUU '^UUi j rp I -^iUU -^iUi j (14) 

\2^010 Xoii J \Xuo XiiiJ 

then the reduced density matrix of |^) is (up to a permutation) Wi © W2 © with 

yyi — In \ l ~ I Tr TS'T Tr TT^ / 



Ai 








Ae, 


A2 








A5, 


A3 








A4 



W2 = (^q' = SS^ + TT^ (16) 

= fn' ^) = S^S + T^T (17) 



It follows from ( fTSi) that the eigenvalues of SS"^ , which are the same as those of S'^S 
can be written as a,Xi — a with < a < Ai; similarly those of TT^and T^T can be 
written as t, Xq — t with < r < Ae- 

The form of f|T6l) and f|T7|) is suggestive of Weyl's problem with A = SS'^,B = 
TT\ C = W2 in the case of ( fT6i) and adjoints reversed for ( |T7I) . Weyl [6l [H] used the 
max-min principle to find necessary conditions 

ai+bi>Ci, «2 + &l>C2, «1 + &2 > C2 (18) 

(with all three sequences in non-increasing order). For 2x2 matrices satisfying Tr ^4 + 
Ti B = TrC, these are also sufficient. We apply Weyl's inequalites to (fT6|) and (|T71) 
and retain the stronger in each pair to obtain 

(x + T > A2 (19a) 
Ai-a + r > A4 (19b) 
a + Ag-r > A4 (19c) 

Adding together the first two inequalities implies 

2r > A2 + A4 - Ai. (20) 

Combining this with 2A6 > 2t and using ([1]) gives 

2(1- Ai) = 2A6 > A2 + I-A3-A1 (21) 

which is equivalent to ([2]). 
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5 Sufficiency 



To prove sufficiency, it suffices to consider a pre-image of the form 

^ = a[(/)i, 02, 03] + ^[01,04,05] + s[(;/)6, 02,04] +i[06, 05, 03] (22) 

and observe that its ffist order reduced density matrix is diagonal in the basis 0^ with 

|a|2 + |6|2 = Ai |sp+|tT = A6. 

Under the assumption that ([T]) holds, the linear relation between the eigenvalues of 7 
and |ap, |tp can be inverted to yield 

|a|2 = i(A2 + A3 - Ag) |6p = i(Ai - A2 + A4) (23a) 
= i(A2 - A3 + Ag) |t? = i(A6 - A2 + A3). (23b) 

With the ordering convention A^ > A^+i, the expressions for |ap, |6p and |sp are all 
positive; and |tp > is equivalent to ([2]). 

In Section [21 we showed slightly more than that ([T]) holds. We also showed that the 
pre-image can always be written in a form in which only four of the coefficients in (1131) 
are non-zero. However, neither of these forms is equivalent to fl22|) with Afc decreasing. 
The equations for the coefficients in one of those forms could have solutions only when 
a stronger inequality than ([2]) holds. In particular, the form obtained from ([7]) in the 
paragraph before (IE]) has solutions only when Ai + A2 < A4 + 1. 



6 General R = N + 3 with N odd 

It is tempting to try to extend the argument in Section [3] to the general case of 
i? = -|- 3 when is odd. Using (jlj) we can conclude as before that 72 must be 
A^-representable with R = N + 2 and thus has an eigenvector l^fi) with eigenvalue 1. 
We can write its pre-image as 

1^2) = am[gi,g2,g3, ■ ..gN-i,gN] + ■■■ + ak[gi, 92,93, ■ ■■g2k-ig2k+2 ■ ■■9n~i,9n] 

+ ... + ai[gi,gi,g4, . ■■gN+i,9N+2] (24) 

where m = j{N + 1) and ak is the coffecient of the Slater determinant which does not 
contain g2k or g2k+i- However, it is not evident that the strong orthogonality condition 
(5^1, $1)1 = holds as was the case for = 3. If we knew that 

\i + {9in9i)<h (25) 

strong orthogonahty would follow, and we could again conclude that 91 is an eigen- 
vector of 7 with eigenvalue {gi, ^ gi) = 1 — \i. 
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For the case A^ = 5,-R = A^ + 3 with odd, Altunbulak and Klyachko [T] have 
shown that Xi + Xr < 1- This is not equivalent to fl25l) because we don't know that gi 
is the eigenfunction for Ar. A condition of the form Xj + Xj' < 1 is sometimes called 
a "strong Pauli condition". When the largest eigenvalue is non-degenerate, we can 
show that strong orthogonality implies a strong Pauli condition with equality. This 
suggests the following 

Conjecture 4. When N is odd and R = N + 3, a necessary condition for pure state 
N-representability of a one-particle density matrix is Xi + Xr = 1, where we have 
assumed that the eigenvalues are in non-increasing order. 

Proposition 5. Let R = N -\- 3 with N odd and consider the decomposition (Hj) o/ a 
one-particle density matrix 7 under the assumption that Ai is the largest eigenvalue. 
Then |$2) has an eigenvector \gi) with eigenvalue 1. If {gi, ^1)1 = 0, then \gi) is an 
eigenvector of with eigenvalue 1 — Ai and this is the smallest eigenvalue 0/7. 

Proof: Let \(pk) denote an eigenvector of 7 orthogonal to both and \gi), and write 

|$i) = a^|0fe®xi) + Vl-a2|^i) 

|'^'2) = bA\gi (S)(f)k^ X2) + Vl-b^ \gi ® ^/'2). 

where we have absorbed any phases into i/jj. Then Xk = Aia^ + (1 — Xi)}P' . Since 
each is strongly orthogonal to \gi) and each is an (A^ — l)-particle 
function with one-rank at most A^. It is well-known [3[ O [10] that this implies that 
iV'j) is a single Slater determinant. Since both l-^i) and \^p2) have one-particle density 
matrices in the same A^-dimensional subspace, it follows that the ranges of these one- 
particle density matrices have a non-zero intersection. Let |/) be in this intersection. 
Then 

(/, 7/) = (1 - a')Ai + (1 - - Ai) = 1 - Afc (26) 

Thus, if, Afc < 1 — Ai, then (/, 7/) > Ai contradicting the assumption that Ai is the 
largest eigenvalue. QED 

7 Further connections with Weyl's problem 

Now assume that gi is strongly orthogonal to $1 and, as in (JTj), write 

7 = Ai|0i)(0i| + (1 - Xi)\gi){gi\ + Ai7i + (1 - Ai)72. (27) 

The A^-representability problem in this situation is reduced to finding conditions which 
ensure that a density matrix is a convex combination of two (A^ — l)-representable 
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density matrices of rank + 1 which satisfy an additional orthogonahty constraint. 
Write 

I'^'i) = ^ Xk^k2...k^_^[gi, g2, ■ ■ ■ gw-i] 

fci<A:2<---fc]v-i 

|'^'2) = 5Z ykik2...kN-i-[9l,92, ■ ■ - gN-l] 

A:i<A;2<...fc]v-i 

Let X, Y be the corresponding anti-symmetric tensors, and let 

XY^ = Xk.M.-kMVk^M.-kM (28) 

k2M---kM 

denote contraction over ■ ■ kM- Then, we can rewrite (I7j) as 

7 - XMicpil - (1 - Xi)\gi){gi\ = xx^ + yyt (29) 

with the constraint ($i,$2) = Tt XY^ = 0. This is a constrained version of Weyl's 
problem. If the R = N + 3 problem could be solved in this way, then by particle-hole 
duality, we would also have the solution to the 3-representability problem. Although 
we we do not know if strongly orthogonality of gi to $i holds in general, this viewpoint 
provides a connection to Weyl's problem that is more general the situation for which 
it was used in Section HI 

For general R (or for R = N + 3 without the simplification that leads to ([7])), 
Coleman's Lemma [1] gives a constrained version of Weyl's problem with 71 = XX"^ 
and 72 = YYK But now 71 is (A^— l)-representable and 72 is A^-representable and the 
orthogonality condition must be translated to tensors of different size. Neverthe- 
less, it now seems clear that what Coleman referred to as a double induction lemma, 
was a constrained version of Weyl's problem. The solution to Weyl's problem was 
given less than 10 years ago, with more recent refinements [9]. Thus, it is not sur- 
prising that the pure state A^-representability problem also resisted solution and that 
Klyachko succeeded by using powerful techniques associated with Schubert calculus 
to solve both problems. 

Acknowledgment: It is a pleasure to recall that most of this work was the result of 
discussions with R.E. Borland and K. Dennis during a visit to the National Physical 
Laboratory in Great Britain in the fall of 1970. 
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